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ABSTRACT 

We use the WMAP maximum entropy method foreground emission map combined with previously determined 
distances to giant H II regions to measure the free-free flux at Earth and the free-free luminosity of the galaxy. 
We find a total flux f^, = 54211 Jy and a flux from 88 sources of fi, = 36043 Jy. The bulk of the sources are at 
least marginally resolved, with mean radii ^ 60pc, electron density rig ^ 9cm"^, and filling factor (pnii ~ 0.005 
(over the Galactic gas disk). The total dust-corrected ionizing photon luminosity is 2 = 3.2 x 10^^ photons s"', in 
good agreement with previous estimates. We use GLIMPSE and MSX 8 /im images to show that the bulk of the 
free-free luminosity is associated with bubbles having radii r 5 - lOOpc, with a mean ^ 20pc. These bubbles are 
leaky, so that ionizing photons from inside the bubble excite free-free emission beyond the bubble walls, producing 
WMAP sources that are larger than the 8 /im bubbl es. We sug gest that the WMAP sources are the counterparts of 
the extended low density H II regions described by Mezged ([19781) . Half the ionizing luminosity from the sources is 
emitted by the nine most luminous objects, while the seventeen most luminous emit half the total Galactic ionizing 
flux. These 17 sources have 4 X 10^' S-' < g < 1.6 X 10^2^-1^ corresponding to 6 x lO'^M© < M, < 2 x IO^Mq; 
half to two thirds of this will be in the central massive star cluster. We convert the measurement of Q to a Galactic 
star formation rate M» = 1 3Mq yr"' , but point out that this is highly dependent on the exponent F « 1 .35 of the 
high mass end of the stellar initial mass function. We also determine a star formation rate of O.14M0 yr"' for the 
Large Magellanic Cloud and O.O15M0 yr"' for the Small Magellanic Cloud. 

Subject headings: Galaxy: fundamental parameters — ISM:bubbles — (ISM) H II regions — stars:formation 
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1. INTRODUCTION 

The star formation rate (SFR) of the Milky Way Galaxy is 
a fundamental parameter in models of the interstellar medium 
and of Galaxy evolution. The rates at which energy and mo- 
mentum are supplied by massive stars, which are proportional 
to the star formation rate, are the dominant elements driving the 
evolution of the interstellar medium (ISM). The hot gas compo- 
nent of the ISM is contributed almost exclusively, in the form of 
shocked stellar winds and supernovae, by massive stars, whose 
numbers are also proportional to the star formation rate. Fi- 
nally, the amount of gas in the ISM is reduced by star forma- 
tion, as the latter locks up material in stars and eventually in 
stellar remnants. Since the star formation rate is of order a so- 
lar mass per year, and the gas mass is roug hly IO'^Mq, either 
the gas will be depleted in 10'' yr, or it will be replaced from 
satellite galaxies or the halo surrounding the Milky Way. 

Estimates of the SFR generally rely on measuring quantities 
sensitive to the numbers of massive stars, including recombina- 
tion line emission (Ha, [Nil]), far infrared emission from dust 
(heat ed primarily by ma ssive stars), and radio fre e-free emis- 
sion. iMezgej (119781) and lGusten & Mezgej (Il982l) showed that 
the latter is dominated not by classical radio giant H II regions, 
but rather by what Mezger called "extended low density (ELD)" 
H II emission. In fact, only 10-20% of the free-free emis- 
sion comes from classical H II regions — the bulk comes from 
the ELD. Free-free emission from H II regions or the ELD is 
powered by the absorption of ionizing radiation (photons with 



energies beyond the Lyman edge, i.e., greater than 13.6eV). 
Thus the free-free emission is often characterized by the rate Q, 
the number of ionizing photons per second needed to power the 
emission (the conversion from free-free luminosity L^^ to Q is 
given by eqn. |7]below). Previous measurements of Q are given 
in table 1, along with the value determined in this work. The 
average of the previous values is Q = 3.2 x lO^-'s"'. 

The io nizing flux can be esti mated from recombination lines 
as well. iBennett et al.l d 19941) use obser vations of the [Nil] 
205 nm line and find Q = 3.5 x lO^^s"'; iMcKee & Williams I 
(1997) use the same observations to estimate 2 = 2.6 x 10 s ' . 

The nature of the ELD is uncertain; it may be associated with 
H II region s, in which case it is also referred to as ex tended H II 
envelopes dLockmanl 1 9761: IXnantharamaiahl 1 985allbl) . The lat- 
ter author lists the properties of the ELD, based on the emis- 
sion seen in the Y\212a line; for Galactic longitudes / < 40° 
the line is seen in every direction (in the Galactic plane) ir- 
respective of whether there was a H II region, a supernova 
remnant, or no point source. The electron densities are in 
the range 0.5cm"-' < n < 6cm"-'; emission measures were in 
the range 500- 3000pccm"^, with corresponding path lengths 
50-200pc; the filling factor is ^ 0.005, and the velocities of 
the H II regions, when present agree well with that of the H272a 
line velocity. 

We note that iTav lor & Cord ei ( Il993h model the free elec- 
tron distribution of the inner Galaxy with two components, one 
with a mean electron density = 0.1 cm"-' and a scale height 
of 150pc, and a second, associated with spiral arms, having 
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{tie) = 0.08cm"-' and a scale height of 300 pc; both components 
are reminiscent of the ELD. 

We present evidence that the bulk of the ELD is associ- 
ated with photons emitted from massive clusters not previ- 
ously identified. We are motivated by the distribution of free- 
free emission in the WMAP free-free map, shown in figure 
H and by comparison of higher resolut ion radio images, e.g., 
Whiteoak et al. ( 1994); Cohen & Green (2001) with GLIMPSE 
jBenjamin et aL2003 ) and MSX CPrice et al.. 2001) data. 

In this paper, we determine the star formation rate in our 
galaxy using the free-free flux measured by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP). We describe our data pro- 
cessing and source identification and extraction methods in ^ 
By comparing to catalogs of H II regions with known distance, 
we estimate the distance to the WMAP sources in ^ The H II 
catalogs are known to be biased against H II regions at large dis- 
tances; we follow Mezger (1978) and Smith et al. (1978) and 
crudely account for this by calculating the luminosity of the 
nearest half of the galaxy, and then doubling the result to find 
the total L^. In ^we examine GLIMPSE images to solidify 
our identifications; in this process we identify 5-75pc bubbles 
associated with the bulk (> 75%) of the emission. We show 
that the bubbles and the free-free emission are both powered by 
massive central star clusters. We derive the ionizing flux Q and 
the star formation rate of the Galaxy in ^ Half the star forma- 
tion occurs in the nine most massive clusters and their retinue; 
the central clusters have M, = 4 - 10 x lO'^M©. We discuss our 
results in |6] In the appendix we describe the machinery needed 
to convert from ionizing flux Q to star formation rate . 

2. MICROWAVE DATA AND WMAP 

The only wavelength range in which free-free dominates 
the emission from the Galactic plane is in the microwave, be- 
tween 10 and 100 GHz, placing this in the center of the fre- 
quency ran ge of cosmic microwave background (CMB) exper- 
iments ( Dickinson et al.ll2003,) . Synchrotron radiation and vi- 
brational dust emission are also important contributors in this 
frequency range. The free-free emission is characterized by a 
spectral index /3, where the antenna temperature Toe i^'^, and 
/? « 2.1. In contrast, the spectral index for synchrotron radi- 
ation is /3 « 2.7-3.2 and for dust emission 1.5-2. In 
order to isolate the free-free component, some form of multi- 
wavelength fitting technique must be used. 

In order to optimize the WMAP measurements of cosmo- 
logical parameters, the galactic foreground emission had to be 
accurately characterized. This was done using a Maximum En- 
tropy Model, r esulting in maps of th e free-free, synchrotron and 
dust emission (Benne tt et alj|2003al) . 

These models agree with the observed galactic emission to 
within 1% overall, with the individual synchrotron and dust 
emission models matching observations to a few percent. In 
the case of the free-free map, the correlation to the Ha map is 
found to be within 12 percent. This indicates that the MEM 
process i s consistent with Ha where the optical depth is less 
than 0.5 dBennett et al.ll2003al) . 

The WMAP free-free model is the only single dish all-sky 
survey of free-free Galactic emission to date, so it is an attrac- 
tive data base to use to measure the Galactic ionizing photon 
luminosity and subsequently the Galactic star formation rate. 

2.1. Data Processing 



We transformed the WMAP free-free maps from an all-sky 
HEALPix map to multiple tangential projections centered about 
the galactic plane. The antenna temperature was converted into 
flux density using the conversion: 



where v is the frequency of the WMAP band, ks is Boltz- 
man's consta nt, c is the speed of l ight, and A74 is the antenna 
temperature dBennett et al.ll2003bl) . To determine all-sky flux 
statistics, an all-sky Cartesian projection of the free-free maps 
was produced. 

The WMAP beam diameter varies from 0.82 to 0.21 degrees 
from the K band to the W band. As part of the map mak- 
ing pro cess, all bands w ere smoothed to a resolution of 1 de- 
gree ( Bennett et al.ll2003a) . The characteristic size of most H II 
regions is of order the smoothed resolution of the foreground 
maps. Thus we suffer from source confusion from regions with 
small angular separations. We discuss our method of separat- 
ing the confused sources in section IZ21 but argue that in many 
cases, spatially separate H II regions are physically associated. 

2.2. Source Identification & Extraction 

Sources within the free-free maps were identified usin g the 
Source Extractor package from iBertin & ArnoutsI dI996l) . The 
fluxes were measured in the WMAP W band, at 93.5 GHz. 
After an automated search over the entire map, a few sources 
were visually identified and extracted. The measured fluxes are 
isophotal with an assumed background flux level of zero. 

Using this method, the smallest extractable flux is approxi- 
mately 10 Jy, with a number of higher flux objects being unex- 
tractable due to confusion within the Galactic Plane. The small- 
est H II region extracted had a semi-major axis of 0.4 degrees, 
half the ~ 1 ° beam diameter of the WMAP free-free map. In 
total, 88 sources have been identified and extracted. 

We have also use d the two-dimensiona l version of the 
ClumpFind routine by Williams et al.l (1 19941) . finding that the 
sensitivity of the isophote parameter provides unreliably vari- 
able sizes and structures for each of the H II regions. Hence- 
forth, we use the sources found by the Source Extractor 

3. DISTANCE DETERMINATION 

As a first pass at distance determination, we use the source 
list of Russeil (2003), who lists both Giant Molecular Clouds 
and H II regions; only the latter are relevant here. In cases 
where the sources have both a kinematic distance and photo- 
metric distan ce, we u s e the p hotometric distance. 

Table 3 in lRusseill (l2003l) lists 481 H II regions; we find 88 
sources, with a much higher total flux. It follows that we have 
likely confused individual sources in comparison to the lRusseill 
d2003i) list. Thus, we have initially assumed that each of the 88 
sources that we have extracted consists of one or more Russeil 
sources projected onto the same location in the sky. We use the 
following procedure to separate these confused sources. 

First, in 13 cases we have a source where Russeil has none. 
In these cases we inspect either MSX or GLIMPSE images to 
identify likely sources, and use SIMBAD to find any HII re- 
gions at promising locations. For example, we find a source 
at / = 6.3 8°, b = +23°, w ith a flux 246.5 Jy, having no coun- 
terpart in iRusseill (I2003h . We identify this sou rce with the C 
Ophiuchi diffuse cloud, at a distance of 140pc ('Drainel lI986l) . 
and find Q/f = 1 A x 10"*^ s~' from the free-free emission; we 
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use a subscript to denote the origin of the estimated luminosity 
(the conversion from L„ = AnD-fy to Q is given in equation]?]). 
This ionizing photon luminosity is reasonably consistent with 
the estimated stellar rate g* = 1-2 x 10"^^ s"' (Panagia 1973), 
and suggests that ^ 35% of the ionizing photons are absorbed 
by dust grains. 

The most outstanding exa mple of a WM AP source with 
no associated H II region in IRusseill (12003b is that at 81.1°, 
b = 0. 5°. This source was, however, mapped bv IWesterhoutI 
(fT958h . who identified it as part of the Cygnus X region. Ex- 
amination of the MSX image shows that there are two large 
bubbles in the region, one centered roughly on Cygnus OB2, 
one on Cygnus OB9. 

We identify the WMAP source at/ = 81°Z7 = 0.5° with the 
northeastern wall of a large bubble in the Cygnus reg i on. Th e 
bubble contains Cyg OB2 (see also 'Schnei der et al.l (l2006l) ). 
The second bubble lies to the south, and appears to contain Cyg 
OB9. The boundary between the two bubbles is a shared wall, 
which contains R ussei]|(l2003h source 118 at Z = 78.5° b = 0.0°. 
His sources 120 and 121 are in the interior of the northern bub- 
ble, near the center of Cyg OB2. The southeastern rim of the 
southern bubble contains Russell's source 1 15. 

We assign a distance D = 1 .7kpc|H anso nl (l2003h to both bub- 
bles (and to the WMAP sources at / = 76.0, / = 78.6 and 81.1). 
We assign the flux from the WMAP source at I = 76.0° to the 
southern bubble, and that of the source at / = 8 1 . 1 ° to the north- 
ern bubble. The flux from the wall separating the two bub- 
bles we rather arbitrarily split evenly between the two. Split 
this way, Q/f = 1.75 x 10^' s"' for the northern bubble, and 
Qff = 1.04 X 10^' s"' for the southern bub ble. We find a to- 
tal free-free flux in the region of 4033 Jy; IWesterhoutI (fT95l 
finds a total flux of 2520 Jy in "point sources" in the region. 

We argue that the free-free flux from the vicinity of the 
northern bubble can easily be powered by Cyg OB2. Count- 
ing only the O stars with spec t roscop ically determined types 
listed in table five of iHansonl (l2003l) yields 49 O stars with 



2 « 5 X 10^°s"'. More recentlvT lNegueruela et al] (12008') find 



50 O stars, and suggest that there may be as many as 60-70 
in the cluster, allowing for some incompleteness due to the 
strong reddening. This is equal to the number of O stars in 
the Carina region as tabulated by Smith (2006), who also gives 
2* = 10^' s"', which we also adopt for Cyg OB2; the total ion- 
izing flux for the region will be somewhat larger, as there are 
a number of O and Wolf-Rayet stars with projected locations 
inside the bubble but outside Cyg OB2. 

We suggest that there must be a similar number of O stars in 
the interior of the southern bubble as well. 

Retur ning to the distance determinations, if there is a unique 
IRusseill (l2003h source at the location of a WMAP source, we 
use his distance as a first guess; there are 43 such objects, about 
half the sample. As in the previous case, we then inspect ei- 
ther MSX or GLIMPSE images at the location of the Russeil 
source. In some cases we find sources we believe to be better 
candidates than the source in the Russeil catalog 



Finally, in 30 cases, we find multiple [Russeill ( l2003h objects 
in the same direction as our WMAP source. We then assign 
a portion of our measured flux to each of the Russeil objects. 
We divide up the WMAP flux using the excitation parameter of 
each Russeil object. The excitation parameter, U oc fi,D^, com- 
pares the ionizing luminosities of the Russeil objects. Using the 
distances provided by the catalog, we calculate the free-free lu- 
minosity of each Russeil object. The result is a separation of the 



confused WMAP source into individual H II re gions wi t h flux , 
distance, and luminosity corresponding to the IRusseill (|2003|) 
objects. 

Using this method, we are able to assign distances to all but 
2 of the 88 regions. (One of the original 13 missing regions cor- 
responds to the Large Magellanic Cloud; we identified 10 using 
SIMBAD and their distances are given in table 3). We assigned 
the average distance of the known sources to the remaining two 
unidentified sources. 

We list 183 H II regions in table 3. For all confused sources, 
the galactic coordinates, semimajor and semiminor axis sizes 
are for the WMAP source, not the individual H II regions. Maps 
of these regions are presented in figures [T]|2| and [3] The distri- 
bution of free-free luminosities dN/dL of these regions is pre- 
sented in figure |4| 

3.1. WMAP sources, the ELD, and dispersion measures 

The WMAP free-free sources range in radius (or semima- 
jor axis) from 0.4° to 10°. The latter is the fitted radius for 
the nearby H II region S2 64 (around A Orionis) at / = 195.05°, 
b = -\\ .995°, D = 400pc dFich & BUtz|[r984 . A visual inspec- 
tion yields a radius ^ 5° or 35 pc, closer to the radius r = 45 .5 pc 
given by Fich & Blitz. As noted above, the effective beam di- 
ameter for the free-free map is ^ 1°. Six sources have mean an- 
gular radii (the geometric mean of the semimajor and semimi- 
nor axes) smaller than the effective beam radius; these are likely 
to be unresolved. The physical radii range from 6pc for ( 
Ophiuchi to ^ 150pc, with a mean radius (r) ^ 55 pc. We find 
a filling factor (j>Hn « 5 x lO""', where (/)hii is the ratio of the 
(summed) free-free source volume divided by the volume of 
the galactic disk assuming disk radius 7? = 8 kpc and scale height 
H = 200pc. 

The ionizing luminosities fall in the range 10**** s"' < Q < 
1.8 X lQ^^s-\ with (e) = 2 X 10^'. The median g = 2.8 x 
10^" s"' (all these values are uncorrected for dust absorption). 

We can determine the mean electron density for each source 
from the expression 
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4'Kr^a(H+)4>' 



(2) 



where a{H^) = 3.5 7 x 10 '-^cm-^ s ' is the Hydrogen recombi- 
nation coefficient (IOsterbrocklll989l) and 4) is the filling fac- 
tor of ionized gas in a given WMAP source region. The elec- 
tron density ranges fromn^, sa lijr^l'^ cm~^ to n^ « 350"'''^ cm"-', 
with a mean n^ = 90"'/^ cm"-'. The density averaged over the 
disk (i.e., multiplying by the volume filling factor (f)Hn) is 
« 0.05 cm"-'. The typical dispersion measure through a 
WMAP source is DM«500cm"3pc. 

The mean mass of ionized gas in a WMAP source is 3 x 
^^(P^/^Mq; the largest sources, with g « 5 x lO^*"' s"', have an 
ionized gas mass ^ 

The density-weighted scale height of the sources is Hhh = 
145pc. 

Recall that the iTavlor & CordesI (1 19931) model for the inner 
Galaxy had two components, with scale heights of 150pc and 
300pc, similar to the scale height we find for WMAP H II 
sources. The mean density of the WMAP sources, averaged 
over the inner Galaxy (i.e., multiplied by the filling factor <6ffff) 
is ng = 0.05 cm"^, compared to the iTaylor & CordesI (Il993l) 
model values 0.1 cm"^ and 0.08 cm"^ f or the inner annulus 
and spiral arms, respectively. Following iMcKee & Wifliams I 
(il997.) . we identify the ELD (the sum of the WMAP sources) 
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with th e arm and annulus components for the iTavlor & CordesI 
(Il993h model. 

3.2. Accounting for the H II region distance bias, and for 
diffuse emission 

We noted above that catalogs of H II regions are known to be 
biased against distant objects, a result apparent in figure [3] We 
followlMezge^(ll978l) : [Smrth et alJ(ll978l) : lMcKee & Williams I 
(Il997h and account for this by doubling the luminosity of 
sources in our half of the Galaxy. This results in 

L^..,ource. = 1 .2 X lO^^crg s"' Hz"' . (3) 

There is a selection effect against low flux sources (less 
than ~ lOJy), as mentioned above, due to the source extrac- 
tion process. The luminosity of a lOJy source at 15kpc is 
2.6 X 1024 erg s"' Hz"', or Q = 3.5x 10^" s"', about one fifth that 
of the ionizing flux of Carina. Since the number counts of free- 
free sources in ground based surveys do not increase much with 
decreasing flux, such sources do not contribute much to the to- 
tal free-free luminosity of the galaxy. 

On the other hand, there does appear to be a diffuse compo- 
nent to the WMAP free-free sky map (diffuse even compared to 
the ELD). The total flux over the entire sky is f^ = 54211 .6 Jy, 
while that in WMAP sources is 36043. OJy. We give a rough 
accounting of this emission by assuming that it arises from gas 
that has the mean distance of the sources, i.e., we multiply the 
free-free luminosity emitted by the WMAP sources by the ratio 
542 11 .6/36043 « 1 .5 to find our final estimate for the Galactic 
free-free luminosity 

L^= 1.8 X lO^^ergs-^Hz"'. (4) 

4. BUBBLES, H II REGIONS, AND MASSIVE STAR CLUSTERS 

We show in th is section that many of the H II regions listed 
in iRusseill (l2003h and earlier compilations are physically con- 
nected. In particular, when several sources appear within < 1 ° 
on the sky, and have radial velocities within Avv ~ ±10km s~', 
examination of Spitzer band 4 GLIMPSE (8 /im) images reveal 
large (10- lOOpc) bubbles, with the H II regions arrayed around 
the rim of the bubble. We interpret these bubbles as radiation 
and H II gas pressure driven structures powered by a central 
massive cluster. Here we give one example; more will be pre- 
sented in a forthcoming paper. 

4. 1 . WMAP sources are powered by massive star clusters 

There are several arguments that the WMAP sources, and 
their enclosed, apparently empty large bubbles actually contain 
the largest star clusters in the Milky Way. 

The first is the very large ionizing fluxes found using WMAP, 
2 3-10 X 10^' s"', for the top 10 or so sources. These sources 
have WMAP-determined radii of order lOOpc, so either there 
are ^ 3- 10 Carina size clusters all within lOOpc, and dN^/dM 
is very different than we believe, or there is a single dominant 
cluster 

The second argument is provided by the shape of the 
GLIMPSE and MSX 8 /im bubbles inside the WMAP sources. 
The bubbles are elliptical, with axis ratios one to two or so. 
This argues for a single massive cluster, which dominates the 
luminosity of the region. 

The third argument is that many of the bubbles show promi- 
nent pillars pointing back to a single location in the bubble, 
again consistent with a single dominant source. 



Finally, we present a quantitative argument for WMAP 
source G298.3-0.34, showing that there should be a massive 
cluster « 4 x 10'^ Mq providing the bulk of the ionizing ra- 
diation. Along the way we show that the classical giant H II re- 
gions associated with this region are powered by compact star 
clusters with masses Q « 7 x 10^°s"', and « lOjOOOM©. 
The total 2 ss 7.7 x 10^' s"' for the region; we show that this 
most likely arises from a cluster at the location pointed to by 
th e giant piflars in figure|5] near / = 298.66°, b = -0.51°. 

ICohen & GreenI yOOl) have shown that the 8 micron emis- 
sion traces free-free emission well. This allows us to use the 
8 /im images to examine the WMAP sources with much higher 
resolution. 

Figure |5] shows the GLIMPSE image in the direction of the 
WMAP free-free source G298.4-0.4. SIMBAD lists 7 HII re- 
gions within 0.5 degrees of the center of the bubble (at I = 
298.5°, b = -0.556°); we interpret 2MASX J12100188-62500 
to be the same source as [GSL2002] 29, and [WMG70] 298.9- 
00.4 to be the same source as [CH87] 298.868-0.432. The 
five unique sources are marked by circles in figure |5] (see ta- 
ble 2). The H recombination line radial velocities range from 
+16kms~' to +30.3 km s~', with a mean around +23kms~'. 
Given the arrangement of sources around the wall, and the 
range of radial velocities, we interpret the source as an ex- 
panding bubble, with mean rtuhbie ~ 55(D/10kpc)pc and ex- 
pansion velocity ^ 7 km s~'. We interpret the H II regions 
around the rim as triggered star formation. The two largest 
H II regions on the rim, G298. 227-0. 340 and G298.862-0.438, 
have fluxes « 47 Jy and 42 Jy, corresponding to 2 ~ 7.5 x 
10-'*°(D/10kpc)s-i and 6.6 x 10^°(D/10kpc)s-i. The total flux 
from the H II regions on the rim is 1 1 1 Jy, compared to the 
WMAP flux of 313 Jy. We suggest that there is a massive clus- 
ter (Q w 3-5 X 10^'erg s"', orM* « 5 x lO'^M©) in the interior 
of the bubble; the pillars point to the location of the cluster 

We note that even so-called "giant H II regions" are spa- 
tially compact, of order a few to ten parsecs in radius (e.g., 
IConti & Crowthei] (|2004|) ); the two classical giant H II regions, 
[CH87] 298.868-0.432 (G298.9-0.4 here) and [GSL2002] 29 
(G298.2-0.3) are prominent in the 8 fim GLIMPSE image, and 
have r adii 3.8 arcminutes, or ^ 10(D/10kpc)pc in 6cm radio 
maps jConti & CrowtheJl2004l) . Radial profiles from the cen- 
ters of the 8/im sources show the l/R surface brightness pro- 
files expected from point sources; see figure|6] These giant H II 
regions cannot be responsible for the much more extended 8 /xm 
emission seen in figure |5] and plotted in figure |6l Nor can the 
two giant H II regions explain the WMAP free-free emission, 
which has r = 0.9° w 160(D/10kpc)pc for G298. 

The surface brightness profile around the large bubble also 
shows a l/R shape at large radii (r > 0.5°). Inside the bubble 
the surface brightness is generally flat, but with a number of 
peaks, culminating in the large peak at r ~ 0.4°, correspond- 
ing to the bubble wall. The total 8 /im luminosity is domi- 
nated not by the known H II regions, but by the large scale 
emission associated with and surrounding the bubble. Figure 
|6] shows the surface brightness profiles of the two brightest H II 
regions associated with G298; recall that both lie on the rim of 
the large bubble. Both profiles merge into the background at 
r ^ 0.3° w 50pc. The figure also shows the azimuthally aver- 
aged radial surface brightness profile from the putative location 
of the massive cluster at I = 298.66°, = -0.507°. In converting 
from degrees to parsecs, labeled along the top of the figure, we 
have assumed a distance D = lOkpc to the object; D = 1 1 .7kpc 
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for V = +23 km s~' in this direction. 

The figure shows that neither of the classical H II regions can 
be responsible for the large scale 200pc) diffuse emission. 
We say this because the \/R scaling for the smaller sources ex- 
trapolates to a very low surface brightness at > 40pc. It also 
suggests that a much more luminous source must be embedded 
in the bubble interior. The surface brightness of the entire re- 
gion also falls off as 1 /r from the point I = 298.66° b = -0.507°, 
as expected if there is a massive cluster near or at this location. 
It follows that the total 8 fj,m luminosity is at least 3.5 times 
that of G298.9-0.4 (the ratio of the surface brightness at large 
radii in the least squares fits) and 5 times that of G298.2-0.3; 
if the emission associated with the H II regions does not extend 
to the edge of the observed 8 /im emission, their contribution to 
the total flux will be smaller 

The azimuthal averaging leads to an artificially thick bub- 
ble wall; surface brightness measurement along radial lines 
show that the radial thickness of the bubble wall is Ar 
4(D/10kpc)pc, about 20% of the bubble radius. 

We noted above that the WMAP free-free source G298 
has a radius of r w 160(D/10kpc)pc, similar to the radius 
200(Z)/10kpc)pc of the 8/im source we find, once again il- 
lustrating the correlation between 8 /im emission and free-free 
emission. 

The total free-free flux in the region is 312Jy, compared 
to 47.4Jy for G298.2-0.3 (- 1/6 of the total) and 42.4 Jy 
(1/7) for G298.9-0.4; note that these ratios are roughly con- 
sistent with the 8/im flux ratios. We estimate a total flux of 
^ 110 Jy for all the classical H II regions in the area, leav- 
ing 202 Jy, which we attribute to the massive central clus- 
ter. We inferred above that the cluster has an ionizing lu- 
minosity Q = 5 X 10^'(D/10kpc)s~', and a stellar mass M « 
5 X lO'^(D/lOkpc)M0, similar to that of Westerlund 1. 

Thus we have a slightly different interpretation of the ELD 
than iLockman (1976) and Anantharamaiah ( 1985a b), at least 
for our most luminous WMAP sources (recall that these dozen 
or so sources supply the bulk of the ionizing luminosity of 
the Galaxy). In these sources, the majority of the ionizing 
flux is produced by a massive star cluster (M 5 x or 
larger). These clusters excite free-free and 8 /im emission out to 
50-200pc. They have also blown ^ 10- lOOpc bubbles in the 
surrounding ISM, as seen in Spitzer or MSX images. The rims 
of the bubbles contain triggered star formation regions, which 
are younger than the central clusters. Because the triggered 
clusters are younger, and substantially less luminous (typically 
by a factor of five), they have not blown away their natal gas. 
As a result, they appear as very high surface brightness free- 
free sources in classical radio emission catalogs (and as bright 
8 /im sources). 

While these young, compact sources are bright and hence 
easily identified, they are not the source of the ionizing photons 
in the ELD. Instead, the massive central clusters are the source 
of the ionizing photons powering the ELD; ionizing photons 
leak out of the bubbles in all directions, since the bubble walls 
are far from uniform. 

Using this definition of an H II region (treating all the H II 
regions associated with a GLIMPSE bubble as one region) al- 
ters the luminosity function. This new luminosity function is 
shown in the right panel of figure |4] At the high luminosity 
end we find dN/dL ~ /v~''^*"'^, i.e., most of the luminosity 
(and stellar mass) is in massive sources (a = 2 corresponds to 
equal numbers per logarithmic luminosity bin). Half the lumi- 



nosity due to sources is in the 9 most luminous objects, with 
Q > 3.2 X 10^"s"' (not corrected for dust absorption). These 
sources have luminosities similar to that of the Galactic center, 
Li, > 3 X lO^^ergHz"', or g > 7 x 10^' s"'. This corresponds 
to cluster masses Mci > 10^ Mp), ranging up to 2.6 x IO^Mq. 

iKennicuttetal J (119891) survey nearby galaxies and construct 
Ha luminosity functions; they find a range of values for a be- 
tween 1.5 and 2.5, with value s below 2 being slightly more 
pr evalent. .McKee & Wifliams I ( Il997h refit the data presented 
in iKennicutt et alj d 19891) using truncated power law fits, and 
find a lower range, 1.4 < a < 2.3, with a mean a = 1.75 ±0.23. 

5. IONIZING LUMINOSITIES OF H II REGIONS AND THE 
GALACTIC STAR FORMATION RATE 

The emissivity of the free-free flux from an ionizing region 
is given by: 



Jf. 



2Ve*> 



In 



y\-'/^Z\,n,e-'"^/'^gff (5) 



where Z is the charge per ion, T is the electron temperature, 
«(. and n,- are electron and ion density respectively, and gff is 
the Gaunt factor. For a fully ionized H II region, we adopt 
lie = n, and Z = 1. Further, we adopt an electron tempera- 
ture, Te = 7000 K for H II regions, and a Gaunt factor g// = 3.3 
( ISutherlandll998h . At radio frequencies we approximate this as 
e(/ = eon^, where eo = 2.1 x 10"-'''gcm^ s""* Hz~' . 

To keep an isotropic H II region ionized, the total number of 
ionizing photons required is: 



n]a{H^)dV, 



(6) 



where V is the volume of the ionized region. 

The total ionizing luminosity (in photons/s) of a given H II 
region is then 

Q,o, = ^^L, « 1 .33 X lO^^L, s-i . (7) 

Using this expression we find the ionizing luminosity of the 
galaxy, before correction for absorption by dust, is Q,ot = 2.34 x 
10^^ photons s~'. 

The final step is to correct for the effect of absorption by ion- 
izing p hotons by dust grains. Following iMcKee & Williamsl 
(119971) . we multiply by 1 .37, and find 

Q,o, = 3.2x10'*^ photons s"' . (8) 
5.1. Star Formation Rate 



To e stim ate the star formation rate fr om Q, we foUow lMezged 
(fT978l) and lMcKee & WilUams I ([T997h . and use the expression 

where (q) is the ionizing flux per star averaged over the ini- 
tial mass function, and (m») is the mean mass per star, in solar 
units. The quantity (fg) is the ionization-weighted stellar life- 
time, i.e., the time at which the ionizing flux of a star falls to 
half its maximum value, averaged over the IMF; all the aver- 
aged terms are discussed in the appendix. 

All of these averaged quantities depend on the initial mass 
function (IMF) of the stars, in particular on the high mass 
slope r of the IMF, as discu ssed in the ap pendix; as an ex- 
ample, and to fix notation, the ISalpeteii ( Il955i) IMF is given by 
^(m) = mdN/dm = Nm~^ , wifli T = 1 .35.^ 

Using the stellar evolution models of iBressan et alJ (Il993h 
we find (fg) = 3.9 x lO^yr (for T = 1.35). This is slightly 



6 



Murray & Rahman 



longer than the ionizing fl ux-weighted main sequence lifetime 
(f„„) = 3.7 Myrs used by iMcKee & Williams I d 19971) which 
is in t urn somewhat larger than the 3 Myrs used by iMezged 
([T978h . This value is only weakly dependent on F. 

The mean ionizing flux per solar mass, {q) / (m*), is much 
more problematic; it depends sensitively on F. Figure|7^ shows 
{q) I {m) using Q{m) as det ermined by Martins et al. (2005) (the 
solid line) and as given bv lVacca et alj (11996) (their ev olution- 
ary m asses); in making this figure we used the Muench et alJ 
( I2OO2I) IMF. The difference between the two estimates for Q{m) 
results in a difference in {q) / {m) of ~ 10%. The filled square 
represents our favored value 



('«*) 



:6.3 X 10'**'s 



■'at' 



(10) 



at r= 1.35^ 

For flie iMuench et all (l2002h IMF (m*) = 0.71 when F = 
1 .35; {q) =4.5 X 10^^ s~'. T his is a factor 5 larger than the value 
quoted bv lMcKee & WilUa ms ( 1997), {q) = 8.9 x 10"*^ s"' ; this 
difference is not primarily a result of our us ing different expres- 
sions for Q{m), since the dashed line uses IVacca et al.l(ll996h . 
as did lMcKee & Williams I (Il997h used. 

We show that this factor of 5 arises mostly from the use 
of a different IMF, with two contributing factors, the use of 
a different value of F, an d a different IMF shape, so that 
lMcKee& Williams I (1 1997h finds fewer massive stars at a fixed 
value of m, even when F is chosen to be the same for the two 
IMFs; in this comparison, we choose F = 1.5 to match their 
work. 

Figure [TJ' shows the mea n ionizing flux per solar ma ss for 
the Scalo-type I MF used b y'McKee & Willi ams I d 19971) (dot- 
das hed Une), the iMuench et al. (2002) IMF (sohd line) , and 
the IChabrieJ (l2005l) IMF (long-dash line), all as a function of 
F. In making this plot we hav e used the relation between Q 
and evolutionary mass given bv lVacca et al. ( 1996h. so that the 
dot-dashed curve goes through the iMcKee & Wilhams 1 (fT997l) 
result. 

From this plot we can see that the variation in F is responsi- 
ble for about a factor of 2 out of the total factor 5 difference; the 
rest comes fr om the different shap e o f the IMF , with the more 
recent IMFs (IMuench et alJ (l200l or lChabried (HOO?)) having 
many fewer low mass stars, or alternately, more high mass stars, 
even for fixed F. 

The figure shows that small changes in F lead to large 
changes in the inferred star formation rate. Recent observations 
of young massive clusters have suggested th at F varie s from 
the Salpeter value (Stolte et al. 2002; Haravama et alj,200 8); 
if confirmed, these variations, combined with the results pre- 
sented here, would lead to large variations in the estimated star 
formation rate of the Galaxy. 

Using the ionizing flux given by Martins et al.l ( l2005l) . we can 
integrate over a .Muench et al.1 (l2002l) -like IMF ;eqn. with 
F as a free parameter. In the appendix we find 



63 X 10*** (m, 



i.35-r\ 



(11) 



where otq « 35Mq is the location of the break in a powerlaw 
fit to Qim) (figure 0). 

Finally, we find a star formation rate for the Milky Way of 

M* =4.1 X 10"^'*e= 1.3M0yr-'. (12) 
Using the iMcKee & Wilhams I ( Il997h value of F = 1 .5 results 
in M* = 2.2Mq yr"' , lower than their 4.OM0 yr"' due to the dif- 
ferent form of the IMF (aside from the high mass slope F) and 



our use of the iMartins et alj (l2005l) tempera ture scale; as seen 
in figure|7] using their IMF and I Vacca et all ( 1 996 ). we recover 
M* w 4M0yr"'. Using the IMuench et alJ (120021) slope, the re- 
sult is 0.9Moyr"'. 

5.2. The Magellanic Clouds 

We were able to measure the free-free flux of the Large 
Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC), 
and thus can provide a star formation rate for each of these 
galaxies. We find = 92.2 Jy for the LMC and = 6.4Jy 
for the SMC. We adopt a distance to the LMC of D = 
48. 1 kpc dMacri et all l2006h. and D = 60.6kpc for the SMC 
(iHilditch etalJl2005h . This leads to free-free luminosities of 



= 2.54 X lO^^ergs-^Hz"' and = 2.81 x lO^^ergs"' Hz"' 
respectively. Using eqns. |7] and [12] we determine a SFR of 
O.UMgyr-' for the LMC and O.O15M0yr-' for the SMC. Our 
estimate for the LMC is slightly lower than but consistent with 
th e estimate of () . 25Mf7) yr~^ found using Ha and MIPS data 
bv lWhitnev et alj (l2008l) . Our estimate for the LMC is signifi- 
ca ntly lower than the Ha estimate of 0.08 Mq yr"' determined 
bv lKennicutt & Hod ge ( 1986) a nd flie IR estimate of 0.05 Mq 
vr~' determined bv lWilke et al] (|2004 . 



6. SUMMARY & DISCUSSION 

We have combined the WMAP free-free map with previous 
determinations of distances to H II regions to measure the ioniz- 
ing flux of the Galaxy. We find 2 = 3.2 x 10^^ s"' , in agreement 
with previous determinations. We found 88 sources responsible 
for a flux of 36043 Jy, out of a total flux of 5421 1 .6 Jy. 

The mean WMAP source radius is ^ 60 pc. Inspection of 
Spitzer GLIMPSE images and MSX images shows that dif- 
fuse 8 i-im emission, which closely tracks the free-free emission, 
gives sizes consistent with the WMAP sizes, e.g., figure|6l sug- 
gesting that many of the WMAP sources are in fact resolved. 

The mean source electron density is 9.3cm~^; hence the 
mean dispersion measure across a source is DM « 540cm"-' pc; 
from figure (Tj most of the sources are within ^ 60° of the 
galactic center. Thus we identify these sources with the inner 
galaxy and spiral arm components of the free electron model of 
iTavlor & CordesI ( Il993h . The density weighted scale height of 
the sources is 114pc. The total volume filling factor of of the 
sources is ~ 0.005. Thus the Galactic mean electron density is 
(n,) « 0.045 cm-3. 

We used GLIMPSE and MSX images to study the WMAP 
sources with higher resolution. We found that the bulk of the 
Galactic star formation (of order half) occurs in ^ 10 sources, 
with g ~ 5 X 10^' s"'. The 8^m images revealed large bubbles, 
with r ^ 20pc, ranging up to lOOpc, in most of these sources. 
We showed that classical giant H II regions associated with the 
WMAP sources were located in the bubble walls, and inter- 
preted them as triggered star formation. We argued that the 
bubbles are powered by massive star clusters responsible for the 
bulk of the ionizing flux in each WMAP source. We estimate 
that these clusters have masses w 4 x IO'^Mq or larger 

We note that there are now a number of slightly older (but 
still young, 10-20 Myr old) Milky Way clusters known to 
have masses in this ran ge; examples inc lude Westerlund 1 
with M 5 X lO'^M© (Brandneret al. 2008), the Arches clus- 
ter M ^ 4 x 10*Mq (Figer et al. 1999) , and the red supergiant 
cluste rs RSGCl nearG25.25-0.15M~ 3 x lO'^Mg dFiger et all 
'200 ^, RSGC2 (/ = 26.2° b = -0.06°) with M - 4 x IO^Mq 
(.Davies et all l2007l) . and RSCG3 (/ = 29.2° b = -0.2°) with 
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M « 3 X IO'^Mq dClark et al.ll2009h . In a forthcoming paper we 
show that almost all of our high luminosity WMAP sources, as 
well as the less luminous sources, are associated with large bub- 
bles seen in GLIMPSE images, and most have faMy compact 
cl usters in the bubb le int erior. 

iLockmanl (1 19761) and lAnantharamaiahl (Il985allbh suggested 
that the ELD, which we identify with the WMAP sources, and 
which accounts for the bulk of the free-free emission in the 
Galaxy, arises from ionizing photons that leak out of H II re- 
gions. We agree that the ELD is closely associated with giant 
H II regions. However, as figure|6]shows, the bulk of the ioniz- 
ing flux powering the ELD arises from massive clusters in the 
centers of large bubbles; the giant H II regions are due to smaller 
(but still large) clusters located in the bubble walls. The mas- 
sive clusters are not readily identified in free-free maps because 
they have blown away their natal gas, and so do not produce 
any high surface brightness emission (either free-free, 8 /^m, or 
even far infrared). 

Using recent estimates of Q{m) and the initial mass func- 
tion ^(m) of stars, we found a Galactic star formation rate of 
M^, = L8M0yr"'. This is somewhat smaller than past deter- 
mination of the Galactic star formation rate. We showed that 



all estimates based on measurements of ionizing radiation are 
highly sensitive to the slope F, where ^(m) ~ at high mass. 
In this case, "high mass" is the critical mass mg « 40Mq where 
stellar luminosities approach the Eddington luminosity. Our 
quoted value of assumes that F = 1 .35, the Salpeter value. 
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APPENDIX 

INITIAL MASS FUNCTIONS, IONIZING FLUXES AND IONIZING LIFETIMES 

We collect here the machinery needed to calcula te the star formation rate from observations of free-free radio emission, following 
ISmifli et al.l (1 1 9781) and iMcKee & Williams I (Il997h . 

We use four initial mass fun ctions, all written in terms of stellar mass in units of the Solar mass, m = M/Mq, with 0.1 < m < 120. 
The first is the lSalpeted(ll955h IMF. 

^(m) = mdN/dm = 7V(F)m"^. (Al) 

Salpeter found F = 1.35. 

The second IMF is the iMcKee & WiUiams I ([T997h version of lScalol([T986l) . which at the high mass end looks like the Salpeter 
IMF, 

mdN/dm = Nm~^ , (A2) 

with TV = 0.063Cf; they ta ke Cf = 1 .4. iMcKee & WiUiams I (Il997l) use F = 1 .5. 
Third, we use a modified lMuench et all (l2002h IMF: 

l' niu > m> nil 

^M,m, (m) = m— = No I mi > m > mz (A3) 

dm ' • — — 



(0.15-r) (-0.73-0.15) 73 ^ ^ 

m\ 'm m2> m > mi 



iMuench et al.l (l2002h found F = 1.21 for the Orion region. As indicated above, mu = 120 and mi = 0.1. We use m\ = 0.6 as the 
characteristic break mass. 



FinaUy, we use flie lChabriei<(l2005h IMF 



iim)-NA ^-^{-''^^] (A4) 
I 0.446m"'" \ <m<Mu 



We use the normalization 



dm 

am)— = L (A5) 



In that case the ionizing flux per Solar mass is 



where 



(a ) r"" dm / r'" 

- ' Q(m)^(m) / / ^(m)dm, (A6) 



(m*) = / £_(m)dm (A7) 

J mi 



is the mean mass per star. 
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We use both the IVacca et alj d 19961) and iMartins et alj (|2005[) compilations of ionizing fluxes as a function of stellar mass; the 
ionizing flux Q{m) is given per star by both. Since many clusters harbor star s with mass in excess of IOOMq, but neither paper 
models stars with M > SSMq, we have added the result of Marti ns et al. I (l2008h . who find Q « 10^" for each of four WN7-8h stars 
with logL/L© > 6.3 (all of which they model by stars with M > 12QMp,; see their table 2 and figure 2 ). These stars are slightly 
evolved, but stifl very young. Figure [8] shows Q{M^) for both lVacca et aLl d 19961) and lMartins et al.l (l2005h . 

The function Q{m) ^ m* for 15 < m < niQ (where niQ w 40), but Q{m) ~ m'-^ for m > niQ. The integral < 2* > (m) ^ mQ(m)(f>{m) ^ 
^1.65 ^Qj. ^ ^jj^ ^ m° for larger m, indicating that the bulk of the ionizing flux occurs for stars with mass around mg, for all 
of our IMFs. Doing the integrals on the right hand side of eqn. iA6i from nii to niQ gives 

-.<^'>, (A8) 



{m 

which fits the numerical result rather well; it is shown as the dotted line in figure |2l 
The ionizing flux-weighted lifetime of a cluster is given by 



/mi I Jmi 

where t(m) is the main sequence lifetime of a star of mass m. 



r" dm / f " dm 

{tQ) = / Q(m)t{m)^{m)— / / Q{m)i(m)—, (A9) 
Jm, m / Jm, m 
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Fig. 1 . — The WMAP free-free map. Note the ~ 60 roughly spherical sources, which we associate with massive star clusters. 
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The WMAP free-free map showing the sources found by Sextractor. 
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Fig. 3. — A map of the location of the WMAP H II regions using the distances derived from lRusseill j20O3ft as described in the text. The shade of the points 
represents the free-free luminosity of the region, log(L^). On the left, we plot all WMAP H II regions, on the right, just the regions with log(L^) > 25 (2 lO" ). 
A bias against distant sources is apparent. 
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Fig. 4. — a) The distribution of free-free luminosity of tlie WMAP H II regions witiiin tiie galaxy, with the corresponding ionizing luminosity indicated on the 
top axis. The number of sources at low flux is reduced by confusion. The dashed line indicates the half luminosity Une, where the sum of the luminosity of the 
sources to the right of this line is equal to half the total measured luminosity in the galaxy. The slope on the luminous end is (dN /dLv ~ a = 1.9 ±0.1. b) The 
distribution for our clumped sources. The slope on the luminous end is a = -1.7 ± 0.2. 
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- ■ . i ■ 

Fig. 5 . — The GLIMPSE 8 fim image in the direction of WMAP free-free source G298.4-0.4. The large white ellipse shows the WMAP source found by Source 
Extractor. We find a bubble in the GLIMPSE image, which we approximate with the smaller of the two white ellipses, having semimajor axis a = 1370" and 
semiminor axis b = 892". We have set the intensity and contrast to show the faint bubble outline, resulting in saturated H II regions. Large pillars are evident at 
/ = 298.67° b = -0.75° and / = 298.5°, b = -0.35°. Also shown (by white circles) are the H II regions listed in table 2; the velocities range from +16km s"' to 
+30.3. We interpret this as a bubble expansion velocity of ~ 7km s~^ The distance to the H 11 regions is D fs 11.7kpc. 
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Fig. 6. — The surface brightness (MJy sr ') as a function of radius, starting from the apparent location of the cluster (/ = 298.66°, b = -0.507°, thick sohd line) 
and from the two giant H II regions G298.9-0.4 and G298.2-0.3 (thin sohd lines; the upper curve near r = 0.1° is G298.2-0.3). The straight dotted and dashed lines 
are least squares fits for 0.01° < r< 0.2° for G298. 9-0.4 and G298. 2-0.3; the slopes are -0.9 and -0.98 i.e., /(r) ~ 1/r. Extrapolating to ( =2°, the upper limit for 
the luminosity of G298.2-0.3 is ~ 1/4 that of the region as a whole, while that of G298.9-0.4 is 1/5 of the total; the ratios found from the free-free emission are 
somewhat smaller The fact that the entire region has a luminosity larger than that of the brightest classical H II regions, combined with the presence of a diffuse 
eight micron emission region much larger than either H II region could illuminate, strongly suggests the presence of a much more luminous star cluster in the interior 
of the bubble. The interior of the bubble shows little emission, as the gas and dust have been pushed to the bubble wall. 
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Fig. 7. — The ionizing flux per sol ar mass for a massive cluster, plotted as a function of the slope F of the IMF for large s tellar masses ; mtdN I dm^ ~ . a) 
The solid line is the result of using the'Marti ns etalJ<2005D expression for Q(m), while the dashed line is the result of using the'Vaccaetal.Vl99^ Q(m). The dotted 
line is the approximation given by eqn. IA8I The two open triangles are the results of Smith et al. 1 197^ (r = 1 .35) and McKee & Williams 1 1997) (F = 1 .5); recall 
that the latter used the Vacca et al. Q(m). The difference between the soUd and dashed lines is about a factor of 1.13 at F = 1.35, so the mismatch between either 
curve and the trian gles is not due to the use of different Q{m). b) Here we plot {q) / (m) for different IMFs, using the Vacca et al. 1 1996) Q(in). The dot-dash line 
uses the^McKee & Williams Ml 997) IMF, so it goes through the open triangle. The solid line represents the Muench et al. (2002 ) IMF, while the long-dash line is 
for the|Chabriet (2005) IMF — it is almost indistinguishable from the Muench IMF. We note that the ratio between the flux per solar mass for F = 1.5 and 1.21 is 
8.37/2.31 ~ 3.5, for all three IMFs. By itself, this dependence on the high-mass slope F leads to a large uncertainty in the star formation rate as measured by any 
method that counts ionizing photons, e.g., free-free. Ha, or [Nil] recombination. Variations in the low-mass end of the IMF will only add to this uncertainty. 
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10 50 100 

M/Mq 

Fig. 8. — The i onizing flux Q as a function of stellar mass M. The open squares (joined b y a dashed line) show the results using the evolutionary masses of 
IVacca et"arHl99 f^. while the solid squares (joined by a solid line) show the results using those of lMartins et aljf2005r) ; both have been supplemented by the addition 
of a slightly evolved model for a I2OM0 star taken from Martins et al. 1 2008). The slope below Mq x 40Mq for both models is d\nQ/d\nM ^ 4, while that for 
Mq is 1 .7, indicating that the bulk of the ionizing emission for any of our standard IMFs comes from stars with M ~ Mq. 
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Table 1 

Galactic Ionizing Flux Measurements 



Galactic Ionizing Luminosity Q Reference 
photons s"' 



3.0 X 10'^ 1 

2.7 X lO'' 2 

4.7 X 10'=' 3 

2.6 X lO'' 4 

3.5 X 10'' 5 

2.6 X lO" 4 
3.2 X 10" 6 



Note. (I) jMezieJ 19751 

(2)|Gusten ik. Mezgei^r ragg) (3)ISm!th etalllTO I 
(4)Jv lcKee & Williams fl9971) (5) IBennett et all 
<1994i) (6) This work. 
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Table 2 

H II Regions Within 0.5° of / = 298.5°, b = -0.556° . 

Name Galactic/ Galactic ii RA(J20()0) DEC(J2000) Flux vv Distance Ref 

degrees degrees hh:mm:ss deg:mm:ss Jy km s"^ kpc 

[KC97c] G298.2-00.8 298.1869 -0.7821 12:09:03.7 -63:15:46 2.4 +16 10.9 1 

[GSL2002] 29 298.228 -0.3308 12:10:04.0 -62:49:27 47.4 -1-30.3 12.3 I 

[KC97c] G298.6-00.1 298.5589 -0.1141 12:13:12.6 -62:39:39 2.8 +23 11.7 1 

[WMG70] 298.8-00.3 298.8377 -0.3467 12:15:19.9 -62:55:52 16.0 +25 12.0 2 

[CH87] 298.868-0.432 298.8683 -0.4325 12:15:29.6 -63:01:13 42.4 +25 12.0 1 



N ote. — (l)[C aswell & Ha vnegfTWTl) (2)|WllM)n et al. i 19 73). Fluxes for [GSL2002] 29 and [CH87] 298.868-0.432 are taken 
from lConti & Crowthea |2004|) : all others are from lT^aswell & Havne j (19871) 



Table 3 
Identified H II Regions 



1 


b 


Semi-Major Axis 
(Degrees) 


Semi-Minor Axis 
(Degrees) 


Free-Free Flux 


Di.sttince 


Distance Reference ^ '^ 


Free-Free Luminosity 

(prtr H?"' 1 

[ cig s nz ) 


Associated Region 


6.4 


23.1 


2.5 


2.1 


247 


0.1 


-1 


1.90E-f25 


COph 


6.7 


-0.5 


1.8 


1.2 


11 


12.3 


6 


2.00E-f24 




6.7 


-0.5 


1.8 


1.2 


9 


13.6 


7 


2.10E-f24 




6.7 


-0.5 


1.8 


1.2 


6 


16.2 


8 


2.00E-f24 




6.7 


-0.5 


1.8 


1.2 


618 


1.6 


9 


1 .90E-f24 


M8 


6.7 


-0.5 


1.8 


1.2 


8 


12.8 


10 


1.50E+24 




6.7 


-0.5 


1.8 


1.2 


281 


2.5 


11 


2.10E+24 


W28 


6.7 


-0.5 


1.8 


1.2 


165 


2.7 


12 


1.40E-f24 


M20 


6.7 


-0.5 


1.8 


1.2 


16 


13.5 


14 


3.60E-1-24 




6.7 


-0.5 


1.8 


1.2 


76 


4.8 


15 


2.10E-^24 


W30 


10.4 


-0.3 


0.6 


0.4 


104 


4.3 


17 


2.30E-^24 




10.4 


-0.3 


0.6 


0.4 


66 


14.9 


18 


1.80E-^25 


W31 


10.4 


-0.3 


0.6 


0.4 


85 


5.5 


19 


3.10E-f24 




10.4 


-0.3 


0.6 


0.4 


20 


14.0 


20 


4.80E-f24 




14.7 


-0.5 


1.4 


0.5 


43 


4.4 


30 


1.00E-H24 





^References to distances are given to table 3 of lRusseil (200 3D 

''Sources with negative numbers in Column 5 have distances given by references as follows: D IDraind I1986D . Refer to the text for more details. 



Note. — Table 3 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 



